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Dui'ing the 6th quarter of this project, we reported 
an attCTipt to invert the equation for the magnetic potential 
in two dimensions: 
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Tests with simulated data showed that, for satellite 
orbits approximating that of MAGSAT, the resultant matrices 
were unstable; i.e., extremely small variations in the 
measured data led to wild excursions in the computed 
values for the coefficients in equation (1) . 


Work done during the quarter has focused on the components 
of the magnetic field. The repressions for the components 
are obtained by differentiating equation (1) . 
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and the tangential component is: 
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From equation (2) and (3) , it can be seen that if 
Z and H are known on a circular boundary, r = a = 1 and the 
Fourier coefficients are computed for each component as; 
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then the coefficients for the internal and external 
parts of the potential can be computed from the relations 


En = nA„ - (n+l)C„ 
and 


Fn = nBn - (n+l)Dn 


( 4 ) 




G = nB„ + , 
n n n n 


Note that we have ignored the zero order term, assuming 
it will be subtracted from the measurements as part of the 
reference field. 

A somewhat similar relation is given in Chapman and 
Bartels for separating internal and external field effects 
given the vertical component and the potential. 

For satellite measurements, however, the data is not 
obtained on a circular boundary, but rathe r on an orbit of 
varying altitut’e. In this case, the terms involving a 
(radius of the circular boundary which xs arbitrary) and 
r (altitude of the satellite) must be included. 

A solution for the coexlicient of the internal and 
external terms in the potential can be obtained by invert- 
ing the matrix equation 


y = Mx 


where 


y ~ 1, Hi/ z 2 1 Hz f • • • • J 

are the measured data values of H and Z 

x=|ai,Bi/Ci/Di/. . . . Aj^, ^ 

are the coefficients for the potential 
and the matrix to be inverted are 
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To test this method, a set of random numbers were used 
as coefficients for the internal and external parts of the 
potential. Using these coefficients, values of Z and H were 
computed along a simulated orbit path given by; 


ALT = 68000 + lOOsin (0) km 


512 values wore computed around the circle 

Two different least-squares procedures wore used to 
solve etjuation (4). The first, a recursive least-squares 
solution was obtained by updating the equation. 
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In these equations is the updated symmetric covariance 

matrix is the updated coefficient vector estimate and 

is the row vector of function values for each new observation. 
For the examples included here, P^ = cl, where I is the unity 
matrix and C ^ 10*. Considerable errors can result if the 

'\i 

start-up values of are too small. 

o 

The second set of compux • -.ions use the standard least- 
squares solution obtained from 
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In Table 1, values for the coefficients of the first 6 
harmonics were assumed in generating both the internal and external 
terms. Both the recursive and standard least squares methods 
were able to recover the generating coefficients quite well, but 
the recursive least squares is slightly less accurate. 


In Table 2, while coefficients of only the first harmonic 
were used in generating the external part of the field, the 
least-squares routines solved for coefficients of all 6 harmonics 
in both parts. Again the recovery was fairly accurate and the 
computed value for harmonics 2-G in the external part were relative- 
ly small. It should be noted that these coefficients have thvn 
dimension of gamma-kilometers, and the actual gamma computed from 
their value would be very small. 


The results in Table 3 show that, although 9 harmonics were 
used to generate the internal and external fields, a solution 
which looks for only 6 harmonics is not thrown off a great deal 
by the presence of the higher harmonics. 


Finally, in Table 4, to more closely approximate the real 
case, the internal terms are made very large and the external 


small. Here the recursive least-squares was significantly worse 
than the standard least-squares in accurate recovery of the 
external field coefficients. 

During the remaining term of the contract, we will apply 
this technique to removal of external field variations from 
actual MAGSAT data. 
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